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NASA’s current predictive capabilities using the ray tracing program (RTP) are 
validated using helicopter noise data taken at Eglin Air Force Base in 2007. By including 
refractive propagation effects due to wind and temperature, the ray tracing code is able to 
explain large variations in the data observed during the flight test. 


I. Introduction 

N ASA has developed a ray tracing aircraft acoustics propagation program, RTP 1 * , based on geometric acoustic 
theory 2 3 ’ 4 , in order to account for refractive effects due to wind and temperature gradients. Theoretical 
studies using RTP have demonstrated that even slight atmospheric changes can significantly impact resulting sound 
pressure levels on the ground at long distances 5 . Now, in order to ensure correct physical predictions, NASA will 
validate RTP against actual flight data. Level flight data from a Mil MI-8 helicopter acquired on September 8 and 9, 
2007, during the 2007 Flight Test at Eglin Air Force Base, will be used for validation. 

II. Test Site Conditions 

Test range C-72 at Eglin Air Force Base in Florida was used for this experiment. The test range spans a 
longitudinal distance of five miles, and, with the exception of a couple of sharp valleys, the terrain is flat. An aerial 
view of the test site is shown in Figure 1. 

Testing was performed in the early morning hours. Table 1 below shows the runs used to validate RTP. All 
cases are steady level flights at an altitude of 300 feet and a speed of 150 knots, with the vehicle moving toward the 
microphone. While it does not greatly affect the results it should be noted that, due to safety concerns, the helicopter 
did not fly directly over the microphone, but had an offset of approximately 300 feet. The primary difference 
between these cases is the relative direction between the wind and the flight path directions. Two flight path 
headings were used for this study, and are shown on Figure 1. Flight path headings, by convention, are designated 
as zero degrees when aligned with north and incremented clockwise to three hundred sixty degrees. 

Table 1 also indicates the corresponding atmospheric conditions used in the propagation predictions. The 
atmospheric conditions were obtained from weather balloon data. 


Table 1. Summary of runs used for RTP validation. 


Case 

Number 

Run Date 

Run No 

Run 

Start 

Time 

Flight Path 
Heading 

Atmospheric 
Condition 
(see Figures 2 
and 3) 

Atmospheric Balloon 
Data 

(reference time) 

1 

9/8/2007 

501702 

7:11 AM 

100 deg 

1 

7:06-7:35 AM 

2 

9/8/2007 

501704 

7:45 AM 

20 deg 

2 

7:35-7:54 am 

3 

9/9/2007 

502711 

7:00 AM 

100 deg 

3 

6:42-7:16 AM 

4 

9/9/2007 

502713 

7:38 AM 

20 deg 

4 

7:38-8:14 AM 
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Figure 1. Aerial view of Eglin Air Force Base C-72 Test Site. 


TIT. Meteorological Data 

Two weather balloons were located at the test site, one at the east and one at the west side of the test site, as seen 
in Figure 1. A cursory comparison shows similar profdes between the west and the east weather balloons. This 
indicates that the atmosphere remains stratified over the test site, and that data from one of the balloons is sufficient 
to characterize the test site. Since the east weather balloon was disabled for safety reasons during our flights, only 
data from the west balloon are available. The weather balloon obtained profiles of temperature, relative humidity, 
wind speed, and wind direction up to an altitude of 3000 feet. Data were gathered continuously during the morning 
hours, with complete profiles obtained about every 20 minutes. 

The temperature profiles for the four atmospheric conditions referenced in Table 1 correspond to the four flight 
cases. The atmospheric profiles are plotted in Figures 2 (temperature) and 3 (wind). The temperature profiles for 
atmospheric conditions 1 and 2 (Figure 2b) show a slight positive lapse rate (temperature decreasing with increasing 
altitude) between the ground and the source height of 300 feet. Atmospheric conditions 3 and 4 shown in Figure 2d 
are inverted (temperature increasing with increasing altitude) and almost constant, respectively. 

The winds at the test site are shown in Figure 3. By convention, the wind direction is labeled as northerly if the 
wind is blowing from the north to the south, and is designated as zero degrees. The wind conditions shown are 
relatively mild and generally blowing from the east at approximately 100 degrees. The wind speeds at the source 
height of 300 feet are approximately 20 ft/sec for all flights. 
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Temp (degF) 


Temp (degF) 


(a). Temperature profiles for altitudes up to 3000 
feet for conditions 1 and 2. 


(b). Temperature profiles for altitudes up to 300 feet 
for conditions 1 and 2. 




(c). Temperature profiles for altitudes up to 3000 (d). Temperature profiles for altitudes up to 300 feet 

feet for conditions 3 and 4. for conditions 3 and 4. 

Figure 2. Temperature profiles for the four atmospheric conditions. 



Atmospheric 
Condition 1 

Atmospheric 
Condition 2 



Atmospheric 
Condition 1 

Atmospheric 
Condition 2 


(a). Wind speed profiles for altitudes up to 3000 feet 
for conditions 1 and 2. 


(b). Wind speed profiles for altitudes up to 300 feet 
for conditions 1 and 2. 
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(c). Wind directions for altitudes up to 3000 feet for 
conditions 1 and 2. 



(d). Wind directions for altitudes up to 300 feet for 
conditions 1 and 2. 



Atmospheric 
Condition 3 

Atmospheric 
Condition 4 



Atmospheric 
Condition 3 

Atmospheric 
Condition 4 


(e). Wind speed profiles for altitudes up to 3000 feet 
for conditions 3 and 4. 


(f). Wind speed profiles for altitudes up to 300 feet 
for conditions 3 and 4. 




(g). Wind direction for altitudes up to 3000 feet for (h). Wind direction for altitudes up to 300 feet for 

conditions 3 and 4. conditions 3 and 4. 

Figure 3. Wind speed and directions for the four atmospheric conditions. 
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IV. Acoustic Data. 

The data used to validate the ray tracing code is a subset of the data obtained at Eglin Air Force Base in August 
and September of 2007. A Mil MI-8 helicopter was used as the acoustic source in all validation cases presented 
here. Vehicle position, velocity, and heading were recorded for all flights using an onboard GPS system. The 
vehicle is moving toward the microphone at about 150 kts, or 253 ft/sec. 

Acoustic data were acquired from a single, ground board mounted, microphone deployed near the east end of the 
test site as shown in Figure 1. The microphone signal was digitized using a 25 kHz sampling rate, and a low-pass 
filter at 12.5 kHz. The microphone system has a dynamic range of 96 dB which is sufficient to acquire 
measurements well below the ambient noise floor, which was between 30 and 35 dB at the first rotor harmonic of 
16.0 Hz for the test cases shown. At the frequency and distances used for the validation cases, atmospheric 
absorption is negligible and the propagation loss is dominated by geometric spreading. The sound pressure level of 
the first main rotor harmonic received by the microphone is shown in Figure 4 for the four validation cases of Table 
1. The levels shown here have been adjusted to account for doppler effects. 


Case 1 


Case 2 
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Case 4 




Figure 4. Mil MI-8 measured first main rotor harmonic for each of the four flight cases. 


The acoustic source is assumed constant for the purposes of this effort. This assumption is reasonable since rotor 
noise in steady level forward flight remains relatively constant within +/- 5 degrees from the plane of the rotor, 
beyond which the signal amplitude begins to drop rapidly 5 6 . Since the Mil MI-8 helicopter rotor has a built in pre-tilt 
of 4.5 degrees from a horizontal plane (see Figure 5), the useable data must be within 9.5 degrees from the 
horizontal plane. Using the combined 9.5 degrees from the horizontal plane and an altitude of 300 feet, it is 
assumed that the acoustic source is constant beyond a range of 1792 feet. 
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Figure 5. Directivity angle over which the rotor source noise is assumed constant. 


V. Predictions 

Acoustic levels are predicted using the ray tracing propagation code, RTP, for the four atmospheric conditions 
described above. As stated earlier, the vehicle is moving toward the microphone at approximately the same airspeed 
and altitude. For the flight test presented, the primary difference between these cases is the relative difference 
between the direction of the wind and the flight path. The refractive effects are caused by the temperature and wind 
gradients. It is the intent of this section to describe the refractive effects for forward propagating rays predicted by 
RTP for these conditions. 

RTP predicts sound pressure level (SPL) losses as a function of range, which is defined as the distance along the 
ground from a position directly under the source (or aircraft). For convenience, the SPL losses are set to zero at a 
range of zero. The predicted SPL levels are shown in Figures 6a and 6b. Currently RTP predictions do not include 
reflections and therefore stop when the acoustic ray intersects the ground. Points A, B C and D in Figures 6a and 6b 
correspond to the ground intersection points for the furthest ray that intersects the ground in each case. These rays 
will be called boundary rays for reference purposes. The figures also show the RTP prediction for a non-refracting 
constant-atmosphere case. 



- Case 1 
Case 2 

Constant 

Atmosphere 



Case 3 

Case 4 

Constant 

Atmosphere 


(a). SPL Losses for Case 1 and 2 (b). SPL Losses for Cases 3 and 4. 

Figure 6. SPL loss predictions for the four cases. 
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It is noted that the SPL losses do not differ significantly (2-3 dB) up to the point at which the rays intersect the 
ground. This is as expected, since the propagation losses resemble spherical spreading at these short ranges. While 
RTP does not predict sound levels beyond the boundary ray locations, investigation of the manner in which the rays 
intersect the ground, along with reasonable assumptions on the refraction effects of the atmosphere, allows one to 
deduce the acoustic levels in these distant regions. Therefore, it is possible to compare RTP predictions at this stage 
to ensure the refraction effects are correct before proceeding with reflections and more advanced propagation 
effects. 

Rays are classified as either upward or downward refracting when they intersect the ground. When an upward 
refracting ray grazes the ground a shadow boundary is formed, and diffraction into the area beyond the shadow 
boundary results in a shadow region with a rapidly attenuating acoustic signal. Likewise, it is expected that 
downward refracting rays are reflected by the ground, and will subsequently be refracted back toward the ground by 
wind and temperature gradients yielding multi-path regions. 

It only takes slight changes in atmospheric conditions to dramatically change long range propagation losses. The 
ability to observe the ray characteristics, and how they change within several atmospheric layers, is important when 
trying to classify associated regions as either shadow or multi-path beyond the point where the rays intersect the 
ground. Table 2 summarizes the ray properties for the four cases, including upward or downward refracting, the 
corresponding expected region’s characteristics, and the ray details for both shadow and multi-path regions. 

For clarity, an explanation of the ray details in Table 2 is given. For shadow regions, there is one unique ray that 
will graze the ground for each case, and for obvious reasons it is called the grazing ray. Associated with this grazing 
ray is a critical elevation angle that describes the angle between the grazing ray and a horizontal plane at the source. 
Any ray emitted at an angle less than (shallower than) the critical elevation angle will turn before it reaches the 
ground, and any ray emitted at a larger angle (steeper) will intersect the ground, unless turned by an intermediate 
layer. The range for both the upward and downward refracting cases is the point where the boundary ray intersects 
the ground. And finally, the azimuthal extent of the shadow boundary shown in Table 2 is defined by the azimuthal 
angles in the ground plane (which are referenced with respect to conventional flight path directions) that locate the 
beginning and the end of the shadow region on the ground. For further explanation of these characteristics see 
Heath and McAninch 5 . 


Table 2. Summary of ray properties for the four cases. 


Ray Properties 

Case 1 

Case 2 

Case 3 

Case 4 

Refraction 

(upward/downward) 

upward 

downward 

undetermined 

downward 

Region prediction 
(shadow/multi-path) 

shadow 

multi-path 

undetermined 

multi-path 

Upward Refracting Ray 
Details 



critical elevation 
angle for grazing 
ray 

11.6° 

NA 

3.7° 

NA 

range of shadow 
boundary 

2355 ft 

NA 

4079 ft* 

NA 

azimuthal extent 
of shadow region 
(with respect to 
flight path 
direction) 

10°-200° 
(A 190°) 

30°-220° 
(A 190°) 

80°-170° 
(A 90°) 

40°-200° 
(A 160°) 

Downward Refracting 
Ray Details 



range of boundary 
ray 

NA 

9075 ft ** 

NA 

2516 ft 


The upward refracting ray that has a critical elevation angle of 3.7 degrees turns before it reaches the 
ground. The boundary is defined by a downward refracting ray at an elevation angle of 5.4 degrees. 

Rays below an elevation angle of 0.7 degrees turn before they reach the ground. The boundary ray defined 
here, with a range of 9075 ft, is the ray defined by an elevation angle of 0.7 degrees. 
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Ray paths have been traced from the source to the ground for several elevation angles and are shown in Figures 7 
and 8 for cases 1 through 4. Each plot shows several downward propagating rays evenly spaced at 0.1 degree 
increments from either the horizontal source plane for downward refracting rays, or from the critical elevation angle 
for upward refracting rays. The boundary rays for cases 1 through 4 intersect the ground at points A, B, C and D at 
ranges between 2000 and 10000 feet as shown in Table 2. These points correspond with those shown in Figures 6a 
and 6b above. 



11.7 through 13.7 in increments of 0.1 degrees. The 
critical emission angle is 11.6 degrees and the 
boundary ray grazes the ground at a range of 2355 
feet (point A). 

Figure 7. Traced ray paths for case 1 and case 2. 



(b). Case 2. Traced ray paths at elevation angles of 
0.7 through 1.5 degrees in increments of 0.1 degrees. 
The boundary ray intersects the ground at a range of 
9075 feet (point B). 


Case i prediction: Case 1 conditions indicate the formation of a shadow region. The shadow boundary is well 
defined, at 2355 feet, and the adjacent rays behave similarly, as shown in Figure 7a. 


Case 2 prediction: Case 2 conditions indicate the formation of a multi-path region based on the fact that the rays are 
clearly downward refracting as shown in Figure 7b. The rays remain downward refracting and will produce a multi- 
path region beyond the boundary ray intersection range of 9075 feet. 



3.8 through 5.8 degrees in increments of 0.1 degrees. 0.1 through 2.7 degrees in increments of .1 degrees 

The critical emission angle is 3.7 degrees and the first are shown. The boundary ray intersects the ground at 

ray that intersects the ground is at a range of 4079 a range of 2516 feet (point D). 

feet (point C). 

Figure 8. Traced ray paths for case 3 and case 4. 

Case 3 predictions : Case 3 conditions indicate neither a shadow region, nor a multi-path region. It can be seen in 
Figure 8a that the first fifteen rays are upward refracting, but turn away from the ground in an intermediate 
atmospheric layer between 50 and 75 feet. The adjacent rays are downward refracting and intersect the ground at 
ranges of 4079 feet and less. Since this case is anything but “ideal”, one might expect the downstream region to be 
affected by both the upward and downward refracting rays. 
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Case 4 predictions: Case 4 conditions indicate the formation of a multi-path region. The rays shown in Figure 8b 
are downward refracting and very compact, in fact, there are 27 rays (covering the first 2.8 degrees of elevation). In 
addition, it should be noted that this case lies along a flight path angle of 20 degrees which is near the asymptote 
where downward refracting rays change to upward refracting rays, which occurs at an azimuthal angle of 40 degrees 
as stated in Table 2. 

Both shadow regions and multi-path regions have distinct characteristics. Acoustic levels in the shadow region 
would be expected to drop instantly to zero if it were not for diffracted sound leaking into the region. The acoustic 

pressure in the diffracted region decays exponentially 7 as p cc e" a Ar , where Ar is the distance from the shadow 
boundary, and a is a coefficient based on the temperature and wind gradients, the ground impedance, source 
frequency, and curvature of the grazing ray. This exponential decay becomes a straight line when plotting SPL 
versus range as shown in Figure 9. A coefficient of .0004/ft was used as the diffraction estimate for the two 
diffraction cases (cases 1 and 3). The diffraction is assumed to begin around 2500 feet. 



Range (feet) 


— - - - Pred Constant 

Atmosphere 

— - - Diffraction Estimate 

— - — ■ Multi-Path (2516 feet) 

— - - Multi-Path (9075 feet) 


Figure 9. Extended predictions for shadow and multi-path regions. 


Multi-path regions on the other hand are dominated by multiple reflections from multiple rays and acoustic 
pressure decays as 1 /-Jr , or equivalently an SPL loss of 3 dB per range doubling 7 . This loss is shown in Figure 9 
for the two conditions, one beginning at a range of 2516 feet to correspond with point D (case 4), and a second 
beginning at 9075 feet to correspond with point B (case 2). SPL losses for a non-refracting constant-atmosphere 
case at the same source altitude of 300 feet are also shown. 

VI. Comparison of Measured and Predicted 

Sound pressure level losses for the 16 Hz main rotor first harmonic are shown in Figure 10 for cases 1 and 2, and 
in Figure 1 1 for cases 3 and 4. Comparison of the two figures (which have similar flight paths for cases 1 and 3, and 
cases 2 and 4, see Table 1) shows large differences in observed losses for moderate changes in the weather 
conditions. For example, a SPL difference of 24 dB is seen between the shadow region and the multi-path region 
for cases 1 (flight path direction of 100 degrees) and 2 (flight path direction of 20 degrees) for the 1 st rotor harmonic 
at a range of 10,000 feet. On the contrary, a difference between cases 3 (flight path direction of 100 degrees) and 4 
(flight path direction of 20 degrees) is sporadic and shows only a 3 dB SPL difference at 10,000 feet. These 
differences, at least in part can be explained by comparing the ray details predicted by RTP. 

Figure 10 compares the SPL losses for predictions and measurements as a function of range. Since the absolute 
acoustic source strength of the helicopter is unknown, comparisons between the predicted and measured SPL losses 
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are achieved by normalizing the data at a range of 2500 ft. This is valid if the noise source is assumed to be constant 
for the duration of each flight. The predictions are shown as “Predicted Case 1” and “Predicted Case 2” for the short 
ranges, and as “Predicted Diffraction” and “Predicted Multipath (9075 feet)” for the long ranges. Figure 10 also 
contains the non-refracted prediction labeled “Predicted Constant Atmosphere” and the ambient noise floor. The 
measurements for case 1 and case 2 show definite shadow and multi-path regions resulting from refraction. 



Range (feet) 


Predicted Case 1 

Measurement Case 1 

Predicted Case 2 
Measurement Case 2 

Predicted Constant Atmosphere 

Predicted Diffiaction 

— - - Predicted Multi-Path (9075 It) 
Ambient Noise Floor 


Figure 10. Measured vs Predicted SPL losses for Cases 1 and 2. 



Range (feet) 


Predicted Case 3 

Measurement Case 3 

Predicted Case 4 
Measurement Case 4 

Predicted Constant Atmosphere 

— - — Predicted Diffiaction 

Predicted Multi-Path (2516 feet) 

Ambient Noise Floor 


Figure 11. Measured vs Predicted SPL losses for Cases 3 and 4. 


Case The SPL loss measured for case 1 in Figure 10 has a clear and steep decline as predicted, until the signal 
falls to the ambient noise floor, which confirms the presence of a shadow region. The SPL decay within the shadow 
region also begins near the predicted shadow boundary range of 2355 feet and is consistent with the diffraction 
predictions described earlier. 

Case 2 \ Case 2 measurements are shown in Figure 10. An initial SPL loss is seen between 2000 and 4000 feet. 
This decline is steeper than expected for a multi-path region, but by 5000 feet the levels are consistent with multi- 
path characteristics. It should also be noted that for these particular conditions there is little difference between the 
constant atmosphere and multi-path predictions. 

Figure 11 compares predicted cases 3 and 4 to the measurements. Like Figure 10, SPL loss measurements are 
shown as a function of range. The predictions are labeled as “Predicted Case 3” and “Predicted Case 4” for the short 
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ranges, and as “Predicted Diffraction” and “Predicted Multipath (2516 feet)” for the long ranges. The figure also 
contains the constant atmosphere prediction and the ambient noise floor. The measurements do not follow either the 
ideal shadow or the multipath regions for these cases. 

Case 3i The data for case 3, as seen in Figure 11, show a well defined shadow region in the beginning, but the 
shadow region becomes overpowered by a multi-path region just beyond 10,000 ft. Since the predicted rays can 
neither define a shadow nor a multi-path region, it is not surprising to see a mixture of the two effects. 

Case 4: Case 4 does not act like a multi-path region as originally predicted, but instead has a very steep drop in 
the sound pressure levels before it switches to a multi-path region at about 6000 feet (see Figure 11). A definite 
cause for this behavior is unknown, however it is suspected that either the flight path heading may be located along 
the transition zone between downward and upward refraction, or, possibly that the multipath region occurs beyond 
the ranges shown here. 

The physical effects predicted by RTP are seen in the flight data, with the exception of case 4. The comparisons 
between Eglin test data and predictions, where well defined refractive properties exist, are consistent. The location 
of shadow region formation and the regions of multi-path propagation show reasonable agreement. The rapid 
attenuation observed within the shadow region is accurately predicted using established diffraction theory. 

VII. Summary 

Validation of NASA’s current propagation capabilities using the ray -tracing code, RTP, is an important step 
toward developing better prediction methods. 

Incorporating wind and temperature refractive effects into the propagation capabilities was the first of several 
key steps in predicting the measured sound pressure levels seen during the 2007 Eglin Rotorcraft Test. Numerical 
predictions from RTP are presented and agree with measurements for three of the four flight test cases. RTP is able 
to predict most of the sound pressure level variations seen during the flight tests, and thus, serves to partially 
validate the modeled wind and temperature refractive effects. 
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